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It is shown that the effects of particle statistics entail reduction in the number of orbital degrees-of- 
freedom in non-relativistic 2-particle systems from 6 to 5. The effect of redundancy in the description 
of orbital motion is found to be in correspondence to the multiplicative phase factor (—1)^" which 
accompany two-particle interchange, where s is the spin of one particle. 



o 

(N 
D 



Oh 
I 

^ ■ 



oo 

> 
m 

in 

cn 
o 

o 



X 



PACS numbers; 



I. INTRODUCTION 



The Pauli exclusion principle is basic principle in 
physics, in particular it is usually related to the expla- 
nation of the shell structure of atoms, conductivity of 
metals, stability of matter, description of the properties 
of white dwarfs, and other phenomena that are of exper- 
imental and theoretical interest. 

A paper by Berry and Robbins has shown that the 
Pauli exclusion principle may originate due to non-trivial 
kinematics of the electronic spins Later, it was found 
that this construction is not unique 0, and alterna- 
tive constructions of the exclusion principle are possi- 
ble. However, an earlier paper by Berry and Robbins ^3)] 
has considered a model AI — Q spin system in external 
magnetic field, and derived the phase factor (—1)'^ mul- 
tiplying the wave-function of the model system, when 
the direction of the magnetic field B is reversed. It was 
shown that the configuration space of the spin system is 
equivalent to the configuration space of two identical par- 
ticles as constructed in Ref. [4] . The derivation in Ref. [1] 
shows that connection between spin and statistics can be 
derived, rather than postulated. In different papers 
a unique connection between spin and statistics for spin 
5 = bosons was derived, based on the identification 
of the symmetric points (ri,r2) and (r2,ri) in the con- 
figuration space of the two particles and exploiting the 
continuity of the boson wave-function. In a related paper 
, the connection between spin and statistics is shown to 
follow if the description of the dynamics involves explic- 
itly ant i- commuting Grassniann variables. On the other 
hand it is known that gauge structure appears in sim- 
ple classical and quantum mechanical systems [S] ■ It has 
been shown that long range forces in di-atom systems 
are mediated by monopole-like gauge fields [9]. Further 
elaboration by Jackiw [l^, [HI has shown that symme- 
tries of dynamic systems can be affected in presence of 
monopole gauge fields. It is therefore reasonable to look 
for a connection between spin and statistics in simple 
systems with two constituent particles. 

In Ref. p^ . three-dimensional variational equation 
of motion for N Coulombically interacting electrons was 
derived, which is different from the conventional 2iN- 
dimensional many-body Schrodinger equation. Unless 



otherwise stated, we use atomic units (e — rrie = h = 1). 

The non-relativistic variational Schrodinger equation 
of motion for one active electron in presence of identical 
spectator particle, and in absence of external (nuclear) 
forces is given by 12 1 



2 ri2 



V'(ricri,r2cr2) = 0, 



(1) 



where g = 1/2, ri2 = |ri — r2| is the relative distance 
between the spectator electron placed at r2 and the active 
electron located near ri, ai and (72 are the components 
of the electronic spins s = 1/2 on an arbitrary but fixed 
spatial z— axis and A is unknown Lagrange multiplier. 
The variational two-body fermion amplitude is given by 

V'(ricri,r2cr2) = (vac|*(ricri)^'(r2cr2)|*), (2) 

where ^'(rcr) is an anti-commuting fermion field opera- 
tor, |vac) is the vacuum state of non-interacting fermions 
and l^*) is unknown state- vector of the interacting sys- 
tem of two electrons. Since the fermion field operators 
anti-commute, then wave-functions of electronic states in 
Eq.(I2]) are anti-symmetric, i.e. 



V'(ri(Ti,r2f72) = -V'(l"2CT2,riCri), 



(3) 



The subsidiary condition of Eq.([3]) provides the descrip- 
tion of the dynamics of the spectator electron, which oth- 
erwise is undetermined. That is because, if 'ip{ra,r' a') is 
a solution of Eq.(IT|), then the re-definitions 



■(/'(rif7i,r2Cr2) 



E 

A1A2 



[/,^^\^^Mr2)^(riAi,r2A2), 



(4) 



are also solutions of Eq.([T]), where ^^1^2 '■''2) arbi- 
trary 16-element spinor matrix, which can depend locally 
on the position-vector of the spectator electron, i.e. the 
phase and the amplitude of the wave-function are not 
fixed by the equation of motion alone. The physical solu- 
tions of Eq.ijT]) are fixed by the Pauli's exclusion principle 
of Eq.(I2]), which fixes the amplitude, the phase and the 
inter-relation between the components of the electronic 
spins of the two-electron wave-function, which otherwise 
remain arbitrary and undetermined. However, since the 
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state- vector \^) is defined up to a global phase, the rela- 
tion of Eq. ([2]) can be written in more general way as 



■(/'(rif7i,r2CT2) = e*®'i/;(i"20'2,ri(Ti), 



(5) 



i.e. the two wave-functions are identical up to a global 
phase-factor, independently on the fact that the fermion 
field operators anti-commute. In the particular case, 
when = 7:, Eq.Q is reproduced, i.e. the particle statis- 
tics phase does not have a direct physical meaning, unless 
the system undergoes a cycle in configuration space, such 
that 9 represents the phase difference between initial and 
final state wave-functions. 

In the more general case of iV-Coulombically interact- 
ing electrons in presence of external one-body potential 
U{y) 12], Eq.([T]) together with the subsidiary condition 
of Eq. (jSj , generalizes to one-electron equation of motion 



\vl + U{v) + g J2 - Aj ^(ra, {vkak}) - 0, 

(6) 



fc=i 



together with iV! — 1 symmetry constraints 

i^il, 2, ... TV) = r^pij[P{l),P{2), . . . P{N)] (7) 

for anti-symmetry of the fermion wave-function due to 
the Fermi-Dirac statistics, by (i) = (r^cTi) we have de- 
noted the coordinates of the i-th electron. The quantal 
statistical phase- factors rjp = (—1)"^ accompanying the 
inter-change of fermions are +1 if the permutation of 
the coordinates involves even number of transpositions 
and —1 otherwise, and iV! is the number of elements in 
the symmetric group Sn- The subsidiary conditions for 
anti-symmetry define the structure of the Hilbert space 
of iV-electron wave-functions. The Hamilton equation of 
motion, together with the subsidiary conditions that the 
wave-function has to satisfy, correspond to Dirac's for- 
mulation of the constrained quantum dynamics (l3| . The 
constraint Hamiltonian approach, has not found realiza- 
tion and application in solving non-relativistic problems, 
such as the calculation of the energy levels of the light 
hydrogen, helium and lithium-like atoms and ions. The 
purpose of this paper is compute the properties of the hy- 
drogen and helium iso-electronic sequences, and compare 
these results to the experiment. 



A. Theoretical formulation 

The state of the interacting iV-electron system can be 
obtained by solving a set of one-particle constraint equa- 
tions for equal sharing of the total energy by the particles 



where 



Xam^O, a = l,...N 



(8) 



(9) 



are the operators of the constraints, A is uniform La- 
grange multiplier, Us(r) is the potential energy of the ac- 
tive electron in an external field t/(r) and including the 
repulsive Coulombic field of the spectator electrons 



(10) 



and Tab = |ra — denote the relative distances between 
particles. In addition there are N first-class constraints 
for identity of the particle spins, which are 



s(s + l)|V), a=l,...iV 



(11) 



where s = 1/2 is the spin of one electron, but otherwise 
s can be regarded arbitrary. The particle spin operators 
satisfy canonical commutation relations 

[(Sa)i, (S6)j] = i5ab£i'jk{Sa)k- (12) 

The constraints in Eq.® are all second-class, since 



[Xa.Xb] 



rab 



' ab 



Pab ^ 0, 



(13) 



where Pab = Pa + Pb is the momentum of the center-of- 
mass motion of the electrons in the (a, &)-th pair. The 
constraints are asymptotically first-class, since they de- 
cay with inter-particle distances as . The second-class 
constraint system can be viewed as a result of gauge- 
fixing in an extended first-class constraint system with 
gauge invariance. We consider N = 2, but we also con- 
sider = 1, since in this particular case the present 
approach reduces exactly to the one-particle Schrodinger 
equation. 



B. Motion of one free electron 

By neglecting temporary the spin constraint, the 
Schrodinger equation of motion in momentum represen- 
tation is given by 



(14) 



where the momentum p of the particle is a multiplication 
operator. The Schrodinger equation is invariant under 
local change of the phase of the wave-function 



«/(p) 



since the momentum p does not change 

P 



e»/(p)pe-»/(p) 



(15) 



(16) 



The phase of the free-particle wave-function is therefore 
uncertain at each point in momentum space. Apart from 
the local phase invariance, there is additional type of in- 
variance of the Schrodinger equation under s-wave re- 
definition of the wave-function, i.e. 



iPip) ^ iPip) + c{E)S{p - 



(17) 
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does not change the equation of motion. Therefore free- 
particle states are defined up to an s-wave, which is 
a consequence of the identity (p^ — k'^)5{p — fc) = 0. 
Therefore apart from conventional symmetries of the one- 
particle Hamiltonian, Schrodinger equation exhibits two 
additional gauge symmetries, it is invariant under local 
[/(I) phase transformations and under s-wave redefini- 
tions of the wave-function. The Hilbert space of states is 
a quotient space, 



^phys ^/^s— "v 



(18) 



The states that are invariant under s-wave transforma- 
tion of the wave-function are 

IV^phys) - IV^) 0, m - 0|V^) ? = 0, m = 0), 

(19) 

for some IV') G i.e. (V'physl^', ' = 0,to = 0) =0, and 
hence the free-particle states do not contain s-wave com- 
ponent. Two states are gauge-equivalent, if they differ 
only by an s-wave 



(20) 



To obtain free-particle wave-functions, gauge fixing- 
conditions CIV') =0 must be imposed, such that to select 
a representative in each equivalence class. 

We further can exploit the local phase-uncertainty of 
the wave-function. The stationary group of the momen- 
tum p is the S0{2) = U{1) group, i.e. 



i?p(x)pi?p^(x) =P, 



(21) 



where i? is a rotation operator and x is the rotation angle 
about the wave-vector p. This also reflects the operator 
identity p • 1 = 0, where 1 = r x p is the kinematic angular 
momentum. The eigenfunctions of the rotation operator 
— are phase factors 



iAx\ 



A 



(22) 



where A is some angular momentum quantum number 
and X = x(p). Each vector v in the space Tp tangential 
to the momentum can be expanded 



ve{v)de + t;^(p)5^. 



(23) 



where the angles (0, i^) specify the orientation of the 
wave-vector p in a space-fixed frame. The phase angle x 
can be defined by the equation 



tanx 



ve_ 



(24) 



which locally changes in the interval < x ^ 27r. In view 
of arbitrariness of v, at each point p of the momentum 
space, we have the freedom to change locally the rotation 
angle x(p) ~^ x(p) + Q^(p) without affecting physical con- 
tent. Since phase transformations of the wave- function 



■(/'(p) i/;(p)e 



iAx(p) 



(25) 



do not change the equation of motion, and no gauge- 
fixing conditions are imposed, the cyclic angle x is a 
redundant gauge degree-of-freedom and the momentum 
space of the particle looks locally like S*^ x S*^. The phase 
angle x, however can affect the quantization of the orbital 
angular momentum. The changes of the momentum un- 
der infinitesimal variation p — p -I- dp induce Abelian 
[/(I) background gauge potential one-form a; = dp • x(p) 
over the momentum space 

dp • (A(p)|iVp|A(p)) dp • cot 6*, (26) 

which maps the gradient of the phase of the wave- 
function, i.e. measures the phase differences between 
wave-function values at different points. When x is re- 
defined locally 



|A(p)) ^ |A(p))^ 



i/(p) 



x(p)^x(p)+Vp/(p) (27) 



the induced displacement field x(p) transforms as a 
gauge field as it should, since the one-form a;(v) = v • x 
takes values on vectors in v g Tp, which are not observ- 
able. If the momentum p is displaced continuously along 
a closed curve C, the phase-factor V'(x) = expzAx satis- 
fies the eigen- value equation —id^ip{x) = ^ipix) at each 
point p. When the wave-vector returns to its original 
direction, the wave-function is multiplied by a Berry's 
phase-factor [l3| 



V'(x)expi7(C), (28) 



which generates a shift of the angle x ^ X + 7(C')/A, i.e. 
a gauge transformation. The gauge-field exhibits Dirac 
string singularity along the entire z-axis, and can not be 
defined globally. Singularity-free induced displacement 
fields can be defined on two overlapping local patches as 
e.g. [li,[l3 



A cos 6* - 1 
p sin 



if, Rn ■.O<0 < (tt + e)/2 



A cosg + U ^^;(^_£)/2<0<^ (29) 
p sm t) 

where the displacement field x^ is regular on the north 
hemi-sphere , while x'^ is regular on the southern 
hemi-sphere R^ . Near the equator R^ P| R^ , these fields 
are related by local gauge-transformation 



x"^ — > x"^ — ie' 



(30) 



The induced displacement field is not rotation symmet- 
ric, since under rotation dp = n x p it changes form. 
Coordinate frame rotations are supplemented by local 
redefinitions of the phase angle x(p): such that rotation 
non-symmetric terms be compensated, i.e. the following 
equation is satisfied 



nxx nxp.Vf 



-Vp/(p) 



(31) 
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where /(p) is a compensating phase function. The gauge 
field does not change form under Galilei boost transfor- 
mations with parameter i5v, i.e. 



Gx(G-ip) = X, 



(32) 



for G = / + iSv ■ r. The gauge field leads to rotation 
symmetric eff'ects through its induced displacement field 
strength two-form F = dx = {diXj{p) — djXi{p))dp^Adp^, 
its dual vector Fi = eijkFjk or magnetic-like field is 



F = Vp X x(p) = -^p, 



(33) 



Further, the kinematic angular momentum 1 = r x p of 
the particle is supplemented by nonkinematic correction 



1 — —p X d 



dW 
dn ' 



(34) 



where W = An • p is the generator of local phase 
transformations, corresponding to phase function /(p) — 
An • p — n X p • X in Eg. ipTj) . The angular momentum 
Ap is the angular momentum stored in the gauge field. 
U{1) gauge-invariant Galilei boost generator d is 



d = r + x(p) 



(35) 



and r = iVp is the noninvariant canonical operator of 
the position. The operator of the position d is non- 
commutative, and satisfies the relations 

[d^,dj] = -iAsijk^, [di,pj] = iSij (36) 

that conserve canonical commutation relation between 
position and momentum. For spin-less particles A = 0, 
the coordinates commute. It is important however, that 
the Jacobi identity is not satisfied 

[[di,d2],d3] + [[d2,d3],di] + [[d3,dl],d2] = -4^A5(3)(p) 

(37) 

The gauge invariant representation of Galilei boost trans- 
formations is based on the operator 



G(v) ~ exp(iv • d). 



(38) 



which generates transformation of the wave-function ac- 
cording to 

G(v)?/;(p) = exp(iv • d) exp(— iv • r)tp{p — v). (39) 

The product of the two exponentials is easily evaluated 
and given by a line integral 



exp(ivi -d) exp(— ivi -r) = exp I i dq • x(q^ 

\ Jp—vi 



(40) 



connecting the frame p — Vi to p. Application of second 
Galilei transformation G(v2), shows that 



where n(p;vi,V2) = ^ ~ Aft/s is the Berry's phase 
[3], i.e. the solid angle subtended by triangle A formed 
by the vertices of the wave-vectors p, p — Vi and p — 
vi ~ V2, as seen from the rest frame p = of the par- 
ticle. Therefore the composition law of the Galilei boost 
transformations is in general non-associative, e.g. [l5| . 
The associativity of the Galilei boost transformations is 
expressed by the equation 



[G(vi)G(v2)]G(v3) = G(vi)[G(v2)G(v3)] 



(42) 



Taking into account Eq. (PT|) . the associativity of finite 
Galilei boost transformations is restored, if and only if 
the flux F of the displacement field strength two-form 
through a tetrahedron enclosing the wave-vector p by 
the three Galilei transformations, is quantizied accord- 
ing to A = N/2. Therefore quantization of the helicity 
A with half-integer numbers is a consequence of associa- 
tivity of flnite Galilei boost transformations. Canonical 
commutation relations between the components of the 
U{1) gauge-invariant rotation operator [U^lj] — iSijkh 
hold for A — N/2. The square of the angular momentum 
operator in Eq. ([M)) can be written as 



sm 



d 



d 



dip 



+ - cos 6*) 



A^ 



(43) 



and the rotation operator about the 2;-axis is Iz 



A. Angular momentum eigen-functions are determined 
by the equations 

l^\lmA) ^ l{l + l)\lmA) , lz\lmA) = m\lmA) , (44) 

for I — |A|, |A| -I- 1, . . . and —l<m<l. Wave-functions 
given by sectional Wu-Yang monopole harmonics 



Yi^a{0,(p) = {9,(p\lmA) 



(45) 



or more explicitly, these are given by means of Jacobi 
polynomials Pn°''^\z) as 

YirnAie, if) = NUl - z)-(A+")/2(l + ,y(A-m)/2 ^ 



(— A-m,— A+m) 



{z)e 



i{A-\-7n)ip 



(46) 



where z = cos 9 and Nim are normalization constants. 
Components of the gauge-invariant rotation operator 1 
satisfy canonical commutation relations [h^lj] = 'i-Sijklk, 
and therefore are connected to the Wigner's rotation 
functions by 

Yi^A{0,^)^Dl,,{~^,e,^) = {lA\e-''^''e'''ye'^''\lm). 

(47) 

The sign of A, sign(A) — ±1 distinguishes left-handed 
from right-handed rotations, which commute. The wave- 
function is independent, single particle states labelled 
by four quantum numbers 



G(vi)G(v2) =exp[in(p;vi,V2)]G(vi+V2), (41) 



IV') = \ElmA) 



(48) 
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and E is the kinetic energy of the particle. The states 
with |A| = 0, 1, 2 . . . form representation of the rotation 
group of integer angular momentum. For A = 0, they 
reduce to the conventional spherical harmonics i^m(p)- 
The states corresponding to half-integer angular momen- 
tum |A| = 1/2,3/2,... define spinor representations of 
the rotation group. For instance, a doublet of wave- 
functions corresponding to / = A = 1/2 is given by means 
of half angles 

riii(0,^) = -sin^e^^ ri_ii(0,^)=cos^, (49) 

on the northern hemi-sphere of the momentum space. 
Second doublet of wave- functions with A = 1/2 with sup- 
port on the southern hemi-sphere is obtained by spatial 
inversion 9 — tt — 9,ip — > ip + t:. Second group of left- 
moving helicity eigen-states of A = —1/2 is obtained by 
complex conjugation of wave-functions of right-handed 
particle states. 
// (up to here) 



C. Motion of two free electrons. 

The Schrodinger equation of motion for one free elec- 
tron in presence of identical spectator electron is, 

(p?- A)^(riai,r2a2) = 0, (50) 

where pi — —iVn is the momentum of the active elec- 
tron. The interchange of the particles' position vectors 
and spins (rici) (r2cr2) leads to identical description 
of the motion of the spectator electron 

(pi- A)^/'(r2a2,ri(Ti) = 0. (51) 

Since the interchange of particles changes only the sign 
of the wave-function, the comparison of Ea. (j50p with 
Eq. (ISlT) shows that the kinetic energies of the two parti- 
cles are equal, i.e. 

P?IV')=P2l^)=2A|V) (52) 

the particles move such that to conserve identical their 
de-Brogile wave-lengths AdB — 27r/V2A. In addition par- 
ticles exhibit identical spins s, i.e. 

=sl\^) =s{s + l)\yj). (53) 

The effects of particle interchange do not involve ex- 
changing energy and momentum and can be represented 
by rotations of unit vectors pi and P2 along with rota- 
tions of half-integer spins, since magnitudes of momenta 
Pi ~ P2 = a/2A and spins s are not relevant for the de- 
scription of the effect of particle interchange, i.e. if for 
instance the momenta are simultaneously scaled accord- 
ing to pi — > e^pi and P2 — >■ e*p2, then the constraint 
equation remains unchanged. The equation of motion 



for the active electron is invariant under bi-local phase 
change of the wave-function in momentum space 

V'(Pi,P2)->V'(pi,P2)e'^'P^'P^' (54) 

and is invariant under s-wave transformation 

'0(Pi, P2) -iAlPi, P2) + 6ipi - A:)c(p2), (55) 

where c(p2) is an arbitrary function of the spectator mo- 
mentum. For comparison, the conventional two-particle 
Schrodinger equation 

(p? + p2-2£;)^(pi,P2) = (56) 

is invariant under bi-local phase transformation, and sim- 
ilar s-wave transformation 

'0(Pi,P2) ^ ^/'(pi,P2) +c5{pi - k)S{p2 - k). (57) 

We further could separate the orbital from the spin vari- 
ables, by demanding that the total wave-function be an 
eigen- function of total spin S = Si -I-S2, together with its 
projection M onto a space-fixed unit-vector P, i.e. 

S^\i^s)^S{S + lMs), P • S|^s) = A^l^s) (58) 

and the wave-function is 

V'(riai,r2a2) = V'5(ri, r2)Cfj^,,, . (59) 

The Clebcsh-Gordan coefficient changes under inter- 
change of spins (Ji o a2 as 

c!Z..-(-^r-'c!.'!^s... (60) 

and implies that under interchange of spatial coordinates 

V's(ri,r2) = (-l)^7/;s(r2,ri), (61) 

the wave-function is multiplied by the phase-factor 
(— 1)"^. We further change the individual coordinates to 
collective coordinates for the relative r = ri — r2 and 
center-of-mass motion R = (ri + r2)/2. The momenta, 
which are conjugate to these coordinates are p = — iVr 
and P — — iVr, respectively. The equation of motion 
reads 

Qp2 + ip . P + ip2 _ a) V5(R, r) = 0, (62) 

and the boundary condition of Eq.Q now reads 

V's(R,r) = (-l)^Vs(R,-r), (63) 

The inter-change of spatial coordinates r — — r in 
Eq. (15^ leads to the equation of motion for the spectator 
electron 

Qp2-ip.P + ip2-A)Vs(R,r)=0, (64) 
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Since Ea. (|62|) and Ea. ([64|) are satisfied simultaneously, 
we have single first-class constraint on the dynamics 

P • p|V^5> = 0, (65) 

that particles share the kinetic energy in equal way, 
and therefore can not be distinguished. The presence 
of spectator particle is non-trivial, since it constraints 
the wave-function of the two-particle state. Hamiltoni- 
ans in Eq. (p^ and Eq. ([M|) together with the constraint 
of Eg. (1651) are translation and rotation invariant. We fur- 
ther constraint the wave-function to be an eigen-function 
of the conserved momentum P of the center-of-mass mo- 
tion 

V's(R,r)=e^P-^V'S.p(r), (66) 

and re-write the equation for the relative motion as 

{p'~k^)^Ps,p{^)^0, (67) 

where fc^ = 2A — /A. The relative wave-function is 
subject to constraint for anti-symmetry 

V's,p(r) = (-l)^^s,p(-r). (68) 

By neglecting effects of spin, the solution of Eg. ((64)) can 
be written as 

V;p(r) =e*''-'e^-^^(""\ (69) 

i.e. the phase of the unconstrained Schrodinger's wave- 
function is re-defined locally by the constraint for equal 
sharing of kinetic energy. Since the particles are free, the 
phase- function /p(r) is a linear function of the relative 
coordinate r, i.e. /(r) — Ap • r, where Ap is a constant 
vector. The non-trivial solution for the compensating 
vector is Ap = — (k-P)P. The wave-function is therefore 
given by 

V.p(r) =e''[''--(''-P)('-P)l, (70) 

and the kinetic energy of relative motion is e = k^ — 
(k • P)^, indicating that the component of the relative 
momentum k on the direction of propagation P is re- 
dundant. The kinematic constraint for equal sharing of 
kinetic energy annihilates the wave-function 

(^eV-plr) =£P-P^'p(r) =0 (71) 

where e is an uniform gauge parameter, i.e. under the 
transformation 

V'p(r)-^7Ap(r)+,5,^p(r)«^(r + eP) (72) 

the two-electron wave-function remains unchanged. Fur- 
thermore the uniform parameter e can be "gauged" into 
a function e — e(r). The particle identity constraint is 
a generator of canonical transformations of the variables 
in the dynamic system, the relative coordinate is gauge- 
dependent and changes as 

<5,r|Vp) =-ie[r,P-p]|7/.p) =eP|Vp), (73) 



i.e. acquires a longitudinal correction. The property of 
sign-change of the vector r = ri — r2 inter-connecting 
the particles under interchange is not unique, since the 
physically equivalent position vector r + eP — >■ — r -|- eP 
does not change sign, unless the transformation ri -H- Y2 
is supplemented by reversal of the sign of P. How- 
ever, the momentum of relative motion is unchanged, i.e. 
<5eP|'0p) = Oj the interchange of momenta p — — p 
is gauge-invariant transformation. The relative angular 
momentum 1 = r x p is gauge- dependent 

J,(rxp)|V;p) =ePxp|7/^p) ^0, (74) 

i.e. 1 1 -I- P X p are equivalent as operators and can be 
identified, in the same way L = RxP^L-|-pxP can 
be identified. The total angular momentum J = L + 1 + S 
is gauge-invariant. If the position vector interconnecting 
the particles is resolved as 

r==r_L + r||, (75) 

where r|| = P(P • r) is the projection onto the propa- 
gation wave-vector, while rj^ is the rejection, then the 
gauge invariance of the wave-function is expressed by its 
independence on the projection r|| and the rejection co- 
ordinate r_L is gauge-invariant coordinate. 

To take into account more accurately the effect of spin 
S, the operator of the constraint is resolved 

P • p = P{px sin 8 cos <i>+py sin sin(i>+p2 cos 8), (76) 

where (P, 8, <&) are the spherical coordinates of the prop- 
agation wave-vector. The inter-particle position vector is 
further resolved in a local basis specified by the propa- 
gation wave- vector P as 

r = rpep + reee + r$e$. (77) 

The rotation matrix which gives the change of coordi- 
nates is 

rp\ I sin 8 cos $ sin 8 sin $ cos 8 \ ( x\ 
re = cos 8 cos <& cos 8 sin $ — sin 8 y 
r$ y y — sin $ cos <i> j \ z j 

(78) 

Similarly the momentum of relative motion is expanded 
over this basis 

p = ppep -l-peee +P*e$ (79) 
At each point P, cylindrical coordinates are introduced 



P = yr^+rl, tan(^=— , z = rp (80) 

The subsidiary condition for equal sharing of kinetic en- 
ergy takes simple form 

dz'^p,sip,'^,z) = 0, (81) 

i.e. wave-function is independent on the longitudinal co- 
ordinate z, and is an eigen-function of the operator of the 
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momentum = pp with eigenvalue kz = kp = 0. For 
each fixed propagation wave- vector P, the wave-function 
satisfies the equation 

+ -dp + X^dl + k^ ^^,s{p. = 0, (82) 

and _E = fc^ is the kinetic energy of relative motion. Un- 
der TT-re-definition of the fiber angle (p, the particles in- 
terchange, and their wave changes according to 

V'P,s(p,^) = e^^>P,5(p,^ + ^). (83) 

The solution of Eg. (15^ is separable R{p)(f){(p) in cylindri- 
cal coordinates, and periodic Bloch-type boundary con- 
dition in Eg . ((831) fixes the solution as 

V^P,5(P, V, a, a') = JA(Me*^^Cf (84) 

where A = S'mod(2?i) is the helicity, which is analogue of 
a Bloch quasi-angular momentum, Ja(z) are the Bessel 
functions of integer order and tt is a characteristic angular 
inter-change period. The sign of A determines left or 
right helicity eigen-states. The inter-change of particles 
r — > — r can be expressed by 

V'p(^ + ^,a',a) = e"^i;p{ip,a\a) = (-1)2^^p(^, a, a'). 

(85) 

where we have used the symmetry property of the 
Clebsch-Gordan coefficient in Eg. ([50)1 . Eg. ([55)) is only 
a consequence of the solution of the equations of motion, 
i.e. if we project the state vector on exchanged config- 
urations {(f + TT , a' , altp) = {—l)^^{ip,a,a'\'tp), the mul- 
tiplicative particle-statistics phase-factor appears auto- 
matically. Therefore, the particle inter-change is a gauge 
transformation of the fiber angle (p ^ ip + n, which par- 
tially compensates the effect of rotation of the phase of 
the spin wave- function. The "local" -type quantization 
of A by Bloch type boundary condition, is inappropri- 
ate, since the phases of wave-functions evaluated at dif- 
ferent points can not be compared, which reflects the 
uncertainty of the phase of the momentum space wave- 
function. 

/ / (up to here) 

The effect of particle identity is shown to entail reduc- 
tion in the number of the initial six orbital degrees of 
freedom to five. 



D. Two Coulombically interacting electrons in 
absence of external forces 

We consider the problem for the Coulomb scattering 
of the two particles, when there is no source of external 
forces. Furthermore, the spin constraints are represented 
very approximately by Bloch-type boundary condition. 
Since the Coulombic interaction is invariant under 
inter-change of particles, then nothing principal changes 



as compared to the case of motion of free electrons. The 
Hamiltonian of the active electron is 

hr, ^ -I'^l + gr^2\ (86) 

By inter-changing the coordinates ri O r2, we obtain the 
Hamiltonian for the motion of the spectator particle 

hr, ^ -^Vl + 9r^,\ (87) 

where g = 1/2. Using Eq.Q that the fermion wave- 
function only changes sign upon inter-change of particles, 
we obtain that 

/iriV'(ri,i"2) = /ir2V'(i"i,r2) = AV'(ri,r2), (88) 

i.e. the particles are precisely identical, since Hamiltoni- 
ans /iri and h,.^ exhibit common eigenvalue A. There- 
fore, equations of motion are consistent only if one- 
particle Hamiltonians commute with each other, i.e. 
[h{l) , h{2)]\^) = 0, which leads to a consistency condi- 
tion 

(Fi2 -pi -F21 •P2)|V') =0, (89) 

where 

F12 = ^ = -F21 (90) 
''12 

is the repulsive Coulomb force of interaction between the 
two particles. Eg. (l89l) shows that electrons move such 
that to screen (compensate) the excess Coulombic force 
in the direction of the total momentum P = pi -I- P2. 
This result has simple classical analogue, since the above 
equation reads 

Pi • Pi - P2 ■ P2 = 0, (91) 

and the difference of the kinetic energies of the two par- 
ticles is a constant of motion 

|(P?-P^) = 0. (92) 

When this difference is vanishing, the particles are pre- 
cisely identical, otherwise they can be distinguished triv- 
ially. Furthermore, if the active electron at point ri 
changes its momentum due to the Coulomb force of his 
partner positioned at r2, the spectator particle changes 
its momentum in strictly proportional way, such that the 
excess force in the direction of the motion of the center- 
of-mass P is compensated. The supplementary condition 
for screening in Eq.(f89|l of the mutual excess Coulomb 
forces can be written more simply as 

ri2-(Vi+V2)|^) =0. (93) 

By means of Eq. lfSS]) . we also have that [ft,(l) — /i(2)]|i/') = 
0, which is a constraint for equal sharing of kinetic energy 
by the particles, i.e. 

[p?-p2]|V;)=0. (94) 
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The constraint for screening of the repulsive inter-particle 
Columbic force can be viewed as a gauge-fixing condition 
for the invariance generated by the constraint of equili- 
bration of the de-Brogile wave-lengths — p^. Re- 
introducing the collective coordinates for relative r = 
ri — r2 and center-of-mass motion R = (ri + r2)/2, 
with the corresponding momenta p — — iVr and P = 
— iVr. Further, the center-of-mass motion is uniform 
and we impose explicitly three additional constraints 
for the conservation of the momentum of center-of-mass 
P = {P^,Py,P,) motion 



zVr(R|7A) =P(R|^), 



(95) 



otherwise arbitrary. We impose planar two-dimensional 
boundary condition for scattering states as 

V(+)(r)«^j,-(r) + /5(fc,^)^^^^^^, p^oo (103) 

where ■0™'^ is the incident Bloch wave, which is superim- 
posed on out-going scattered Bloch wave J^^'^^ of ampli- 
tude f{if). The planar Bloch wave- functions of electronic 
states exhibit partial wave-expansion as 



rZisikp)- (104) 



i.e. the two-electron state is an eigen-state characterized 
by the momentum P, i.e. 



^ = e'^-^Vp(r). 



(96) 



Similarly, the scattering amplitude exhibits Bloch repre- 
sentation 

/5(fc, ^) = e^'^ J2 e'"'^/2«+s(fc), (105) 



By separating the center-of-mass motion, the Hamilto- 
nian for the relative motion of the two particles becomes 



Hrel = P' + \r\-\ 



(97) 



subject to the supplementary condition for screening P ■ 
r|'0P,s) = and for equal sharing of kinetic energy 



P • p|V'p,s) - 0, 



(98) 



and satisfies 



fs{k,ip) = {-lffs{k,ip + 7r), 



(106) 



i.e. it is symmetric for scattering in a spin-singlet and 
anti-symmetric otherwise. Instead for the physical wave- 
function 4!^~^\kp), we solve this equation for the regular 
wave- function R{kp), which is subject to the boundary 
condition that 



and P is a unit vector in the direction of propagation of 
the center-of-mass motion. In Cartesian coordinates r = 
irp,re,r^) with z-axis parallel to the propagation wave- 
vector, the pair of supplementary conditions become 



0, 



(99) 



which is a pair of second-class constraints, that show that 
the longitudinal relative coordinate rp is locally redun- 
dant. The reduced Hamiltonian for the planar orbital 
motion of the internal degrees-of-freedom simplifies as 



Hrel =P%+P% 



1 



(100) 



The wave-function is subject to the boundary condition 
^s.p{re,r^) = {-lfi^s.p{-re,-r<i>). (101) 

Introducing the cylindrical coordinates tq = p cos ip and 
r$ — p sin if, the Hamiltonian reads 



(102) 



To comply with scattering state boundary conditions, we 
specify the orbital collision plane to be formed by the in- 
cident wave- vector kj and the wave- vector k/ of the out- 
going scattered wave and therefore P = k^ x k/ speci- 
fies the orientation of the orbital collision plane, which is 



limp-l^ii?|Ai(p) = 1. 



(107) 



and satisfies the equation 



■t^Ra{p) + -^r^^^p^ + -] ^a(p) - 0, 

dp Pdp V P P- ) 

(108) 



where A = S" mod2fi. is the heficity, k = ^2(A- P^/S) 
and the solution depends only on |A|, and we further take 
A > 0, which is equivalent to take n > and consider 
right-handed electronic states. Making the substitution 
-Ra — uk/ \/p in Eq. (|108l) . we obtain the Whittaker's 
equation 



1 7] 1/4 -A^ 



4 z 



ma(^)-0, (109) 



where z = —2ikp and rj = —i/2k. The general solution 
of the Whittaker's equation is 



(110) 



where Aa and Ba are integration constants and 
W±,j^A{iz) are the two linearly independent Whittaker's 
functions of second kind. The unknown integration con- 
stants can be obtained from the boundary condition for 
the regular solution Eq.(1107|) . The asymptotic of the 
Whittaker's functions [iM, when \z\ — >■ cxd. 



(111) 
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gives the asymptotic of the scattering-state wave- 
function 



ua(z) ? 
+ BAe 



-i{kp~\og 2kp/k) 



p — > OO, 



(112) 



as hnear combination of irregular Jost solutions, specified 
by the boundary conditions 



lim eT'['=^-'°s2fcp/fc]jr(±)(fcp) ^ 



(113) 



and describe polar Coulomb waves outgoing from [with 
sign (-1-)] or incoming [with (— ) sign] towards the origin 
p = Q. The Jost solutions and Whittaker's functions are 
identical up to multiplicative A-independent constant, 
more specifically the relation is given by 



T^^Hkp)^e^/^''W±r,A±z), 



(114) 



and these functions can be identified by A-independent 
re-definition of integration constants 

A Ae'^/^fe ^ f(-)(fc)^ B Be'^/^fe ^ _f(+)(fc)^ 

(115) 

The wave-function can be re-written as linear combina- 
tion of the two Jost solutions 

u{p) = --^[f(-)(fc)^(+)(fcp)-f(+)(fc)J^(-)(fcp)], (116) 

where w{k) — W^[J^'^^\ = 2ik is the Wronskian of 

the two Jost solutions, i.e. W[f, g] — fg' ~ f'g and prime 
denotes radial derivative. The integration constants f*^^' 
are the Jost functions, which are given by the Wronskians 



f(±)(fc) = M^[J-(±) 



(117) 



which we evaluate at the origin p = 0. By using the 
asymptotic of the irregular solutions near the origin 



lim z^-'/^W„a{z) = ^^^^^ — 



(118) 



and evaluating Wronskians with the help of the boundary 
condition in Eq. (|107p . we obtain Jost functions as 



f(±'(A:) = e'^/2'^(T2ifc)i 



/2-A_ 



r(2A + 1) 



r(A±i/2fc + l/2) 



, (119) 



where T(z) is the Euler's gamma function. Then the 
asymptotic scattering-state wave-function p — >■ oo is 
given by 



p7r/2fc 



{2k) 



-1/2-A 



r(2A-i- 1; 



|r(A+ 1/2 + V2fc)| 



X sin ( fcp - - log 2fcp - — - - 



(120) 



where i5a are elastic scattering phase shifts 

5a(A:) = argr(A-M/2-t-i/2fc), (121) 



relative to the asymptotic of the non-interacting Bessel's 
function J\{kp), 



2 ( ttK IT 

cos Z 7 

-KZ \ 2 A 



(122) 



The Coulombic S'-matrix is given by the factor of the 
Jost functions 



(123) 



The physical wave- function differs from the regular wave- 
function be a normalization constant determined from 
the Jost function, i.e. 



^i-)(M = iVA(fc)"^(^^) 



/2kp ' 



and therefore 

NK{k) = e 



W2fc(gfc)A r(A + l/2 + V2fc) 



(2A)! 



(124) 



(125) 



As a result, we obtain the partial- wave scattering ampli- 
tudes /a as 



-i-iT /A 

/27rfc 



Jg2i(5A(fc) _ ^ 



(126) 



We next evaluate a differential cross-section for scattering 
in a given line segment in the collision plane as 



(127) 



To evaluate outgoing scattered probability flux through 
a solid angle we vary the wave- vector kj, such that 
the unit vector x k/ normal to the scattering plane 
rotates on angle a about the axis of incidence k^, and a 
generated linear flux is 

dP„ = |/((^)|2(sin(^da). (128) 

Therefore elastic scattering cross-section is given by 



da - dP^dP^ = \f{^)\''dn, 



(129) 



where dfi — sin ipdipda is the solid angle of observation 
in the space-flxed reference frame. The differential cross- 
section is given by 



da da^kf k^) 



dn 



= \fik,^)\' 



(130) 



and exhibits characteristic dependence on the fourth 
power of the amplitude. 
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E. Quasi-classical approximation 

In quasi-classical approximation, the two-electron 
wave-function is a phase-factor 



js 



(131) 



expanding the phase S = So+hSi + . . ., to zero'th order in 
the Planck's constant, we obtain the equation of motion 
for the active electron as 



1 



[Vr,^o(ri,r2)]^ 



g\ri -ral 



A, 



together with the constraint for particle identity 

[Vr,S'o(ri,r2)]2 = [Vr,S'o(ri,r2)]2. (133) 

The quasi-classical Coulmbic action has the form 

5o(r,R) = P-R + ap(ri), (134) 

where the action for the relative motion satisfies the equa- 
tion 



with parameter p = and eccentricity e = A = 
\/l + kP'K^ > 1, and the point of closest approach on the 
trajectory {(p — 0) is pmin = p/(e — 1)- The cross-section 
for elastic scattering can be derived from the conservation 
of the Runge-Lenz vector. In a reference frame, where 
the center-of-mass motion is at rest P = 0, we choose 
the direction of incidence be the negative half of the 
X-axis of the laboratory frame, the asymptotic momen- 
tum of relative motion is = kix, the angular momen- 
tum of relative motion corresponding to this choice is 
A = xpy — ypx = —bk, where b is the impact parameter, 
(132) and By is a unit vector in the plane of the orbit. The 
Laplace-Runge-Lenz vector prior to the collision is given 

by 



1 



(135) 



where l/r± is the screened planar Coulomb potential. 
Re-introducing plane polar coordinates rj^ = {p, ip) to 
describe a collision, the quasi-classical equation reduces 
exactly to the Hamilton- Jacobi equation of the planar 
Kepler problem 



[dpcrf + hd^af + - ^ k'' 
P P 



and has solutions 



a = ±Atp + (jk{p), 



(136) 



(137) 



where A labels the helicity. The quasi-classical approx- 
imation holds if A ^ 1, and therefore we will neglect 
effects of quantization of A. The planar Kepler problem 
exhibits dynamical symmetry (e.g. [ISIill)) due to the 
conservation of the planar Laplace-Runge-Lenz vector, 
which is given by (we use the classical expression, which 
does not involve hermitian symmetrization) 



A = p X 1 + f , 



(138) 



where 1 = r x p is the relative angular momentum, p — 
VrO'(rj^) is the relative quasi-classical momentum and 
A • 1 = 0. The classical trajectories of relative motion of 
the two-electrons can be obtained from 



Y ■ A — pA cos If — K 



(139) 



which leads to the conical section equation, specifying 
unbound hyperbolic Kepler orbits 



P 
P 



- 1 + e cos (/?, 



(140) 



bk^e 



and similarly after the collision 

Aout — Bout + bk^Uout, 



(141) 



(142) 



where Bout = f ^ is a unit vector specifying the outgoing 
direction of the scattered particles, hout = P x Bout is 
a unit-vector in the collision plane, and P is the unit 
vector normal to the plane of the orbit. By projecting 
the conserved Runge-Lenz vector onto the direction of 
incidence, i.e. 



we obtain that 



— 1 = cos X — bk^ sin x, 



(143) 



(144) 



where % is a rotation angle, cosx = B^ ■ Bout- From 
Ea. (|144p we obtain the relation b — b{x) between the 
deflection angle and the impact parameter as 



b = 



1 



k^ tanx/2 

Differentiating with respect to Xj we obtain 



db 
dx 



1 



2fc2sin^ x/2 



(145) 



(146) 



The classical differential scattering cross-section is given 
by da — 2'iTbdb, and by using that dU, — 27rsinxdx, to- 
gether with Eq. p45| and Ea. ([T46)) we obtain the Ruther- 
ford formula for the cross-section 



da _ 1 

'dn ~ 4fc4sin'*x/2' 



(147) 



Taking into account reflection symmetry % — )■ vr — x, i.e. 
that we can not distinguish between forward and back- 
ward scattering when the particles are identical, we ob- 
tain that 



da{-K - x) 

dn 



4fc4cos4x/2 



(148) 
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The total cross-section for elastic scattering is obtained 
by the sum of the two contributions 



dn 4fcHsin4x/2 cos4x/2; ^ ' 

The cross-section for scattering at small and large angles 
is highly divergent. The cross-section can be defined only 
when the center-of-mass motion is at rest. The particles 
are identical and have equal kinetic energies pi — p\ due 
to the statistics, and there is no reference frame where 
only one of the particles is at rest. This also means that 
the problem does not exhibit spherical rotation symme- 
try, instead it exhibits cylindrical rotation symmetry. 



II. CONCLUSION 

We show that in the particular case of systems with two 
particles, that the constraints of particle identity entail 
reduction in the number of internal degrees-of-freedom 
from six to five. The effect of redundancy in the descrip- 
tion of orbital motion in the two-particle gauge system 
is found to be in correspondence with the multiplicative 
phase-factor (—1)^, where S = {0, 1} is the total spin. 
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